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Abstract 

We discuss the resummation of the potentially large 0^(/?oa s ) n ) QCD correc- 
tions to the total cross section of the process e + e~ — * tt near the threshold. In this 
approximation, the cross section factorization into the short- and long-distance 
parts is valid. The short-distance correction is reduced to the production vertex 
renormalization. It amounts to —9.5% and is well under control. The long distance 
corrections are accounted for as the effect of the coupling's running in the QCD 
potential. We argue that the accuracy of present predictions for the cross section 
in the scheme with the running QCD coupling is about 10%. 
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1. The theoretical studies of the process e + e~ — > ti near the threshold are of interest 
for many applications. First of all, the top quark will be studied in the near threshold 
kinematics at the Next Linear Collider, and precise measurements of its mass and width 
are among the most important experimental issues |TJ. From this standpoint, reliable 
estimations of the QCD perturbative effects are mandatory. 

The general approach for the problem of the top quark production near the threshold 
was suggested in JJ, where it was demonstrated that the top quark width plays the role 
of an infrared cutoff which suppresses all nonperturbative long-distance effects. The size 
of the nonperturbative corrections was estimated in || M to amount to less than one 
per cent. Therefore, one could hope to predict the total cross section by perturbative 
calculations with such an accuracy. One-loop QCD corrections were discussed in p, ||, 
|6| and Higgs-induced electroweak corrections, which have the same order as the QCD 
ones, were discussed in [0. Recently, the 0(a 2 s ) QCD perturbative corrections to the ti 
threshold cross section were calculated in || |[J and were found to be comparable in size 
with the first order correction. This result can cast some doubts upon the applicability 
of QCD perturbation theory to the problem under discussion. 

In this note we make an attempt to answer the question: "How can higher order 
perturbative corrections modify the NNLO result?". Surely, an exact answer to this 
question deserves an exact knowledge of those higher order corrections which is far beyond 
our reach at the moment. Nevertheless, there exists a method to estimate a contribution 
of higher orders, which is known to work well at different problems (see for example 



|T0| , p] , [121, [13], [14], [15], [16[ and references therein). This method exploits the large value of 
/?o =11 — I^Yf, the lowest order coefficient from the QCD beta-function, and consists in 
an approximation of the complete result for the (n+ l)th order, by its counterpart of the 
nth order in (3q. The latter is found by the naive non-abelianization (NNA) prescription: 
in the (gauge-invariant) result for a set of Feynman diagrams containing n light quark 
vacuum polarization bubbles, one replaces Np by (— 3/3o/2) n . In other words, taking 
into account only the running of the QCD coupling, one assumes it to be a reasonable 
approximation for all radiative corrections. 

A nice feature of this approach with respect to the threshold cross section analysis is a 
complete factorization of hard and soft corrections, both of which reside on its own scale. 

In what follows, we first discuss the intermediate kinematic region a s < u < 1. We 
resum O(^(f3oa s ) n ^j terms to the first two leading coefficients in the expansion in v. Then, 
we consider the cross section near the threshold, where short -distance effects are factorized 
to the normalization of the cross section and the long-distance effects are absorbed to 
the static potential with the running coupling. We perform a numerical evaluation with 
such the potentials to test the influence of the higher order QCD corrections on the cross 
section. 

2. Let us first consider the kinematic region a, < i; <C 1, where v is the relative 
velocity of quark and antiquark. In this region, we have an extra small parameter, a s /v, 
so that the interaction between the slowly moving quark and antiquark can be treated 
perturbatively. In the lowest order, this interaction is caused by the one-gluon exchange 
between two particles. Corresponding one-loop correction to the pair production ampli- 
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tude, A (v), can be easily obtained using the QED result by Schwinger [|T? 



Here the former term is determined by the low-energy scale, so that all powers of Cfol s /v 
should be resummed in the vicinity of the threshold. The latter term has its origin at the 
relativistic scale. 

As discussed above, we would like to consider the effects of vacuum polarization only. 
The most appropriate way to do that is to substitute l/{k 2 — A 2 ) for 1/k 2 in the gluon 
propagator, and to calculate the amplitude A\(v) as a function of the "gluon mass" A. 
Then the first -order in a s correction to the production amplitude, which is valid to all 
orders in /3 a s , equals |L2], [13], [L4], [15|, [Tl[ 



where 

is the one-loop polarization operator, fi is the normalization point for the coupling con- 
stant, a s = « s (/i), while C is a scheme-dependent constant. In what follows, we use 
C m = -5/3. The result for A x (v) reads @, @: 

Ax(v) = -^arctg- 4 

V z 

2C F a s ( 3vr z 6 - 2z A - 2z 2 - 6 z 4l 2 3] 

H . arch 2 mz + z H — > (5) 



3tt [2^ zv 7 ^! 2 2 

= ^a(^) + ^a- (6) 

Here 2 = A/m, m is the mass of the t-quark. Inserting the difference A\(v) — A (v) into 
(0), we obtain the amplitude in the intermediate kinematic region as a sum of the soft 
and hard terms, 

A(v) = A s {v) + A h . (7) 

3. Now we would like to find the hard correction, A , due to (|5]). Inserting (|5|) into 
the integral (0), one can obtain its expansion in power series over a = /? a s (/i)/(47r): 

A ^_^M Er/) (8) 

71 n=0 

where the lower order coefficients at jj, — m are: 

r = 2, n = ^, r 2 = ^- + r 3 w 57.05, r 4 w 1131. (9) 

o 3 18 

The zeroth order coefficient coincides with that from Schwinger 's result ([!]), the first order 
one is in accord with the result obtained in ||19||, while r 2 , r 3 and r 4 are new. Due to the 
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Figure 1: A% exact »(a) (solid line); A h (dashed and dotted lines) expanded in a n , as a 
junction of a. Vertical lines mark the actual values of a at the threshold for top (m = 175 
GeV), bottom (m =4-8 GeV), and charm (m = 1.4 GeV) quarks. 



fact that the vector current's anomalous dimension is zero, one can easily restore the in- 
dependence of the results. In Fig.l, we compare the results of the lower orders inclusion 
with the exact (in NNA approximation) result. The actual values of a for c-, b- and 
t-quarks are indicated in order to demonstrate how the perturbation theory works for the 
current renormalization at the corresponding thresholds. 

Taking a s (m t ) = 0.1, so that a = 0.061, we can calculate an "exact" value of the hard 
correction to the amplitude of tt production: 

A h w {a) Uo.061 = -2.236^ (10) 

7T 

« (-2-1.33a-41.7a 2 )^^. (11) 

7T 

The second equation is an extrapolation formula which works in the region a = (0,0.1) 
and can be used to obtain A« exSuCt ,, (a) at different a s {m t ). Adding up this result to unity, 
we obtain the renormalization of the vertex 7*tt at the threshold, by the hard QCD 
correction (in NNA approximation). The O{a 2 s (3o) correction amounts to about 6% of 
the one-loop 0(a s ) result Aq, the O^otlfll) correction gives about 4% and the rest in 
the sum gives 2%. The sum of higher order corrections constitutes 12% of Aq. The net 
QCD correction to the renormalization of the 7*tt vertex amounts to —9.5%, in the NNA 
approximation. 

We would like also to compare this number with the two-loop QCD correction to 
the vector current normalization constant, ^wo-ioopt^fact)) derived in [^, ||, ||, which 
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depends on a factorization scale. We take eq.(22) from (or eq.(39) from 0), to obtain 
Awo-ioo P K) = -2-34^ and ^ wo _ loop (m 2 /2) = -2.14^. One can see that our 
result (pT0|) is rather close to these numbers. 

We conclude that the discussed corrections to the short-distance part of the production 
amplitude are well under control. 

As for the long-distance part in the intermediate region, we obtain 

A\v) = °*^'Et n a», (12) 
where the lower order coefficients are: 

1 D D 2 7T 2 

to = -, tl = - T> t 2 = T + T , (13) 
D 3 „_ D A „ 8tt 4 



£3 = 7T 2 D, t A = + 2L>V 
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and D = C — 21n( — ). The zeroth order coefficient coincides with that from Schwinger's 
result (|TJ), the first order one coincides with the result obtained in [O], while t 2 , t 3 and t± 
are new. 

It is worthy to stress that the intermediate region has only academic interest, since at 
a s = 0.1 there is no room for the strong inequality « s < f < 1. Either 0(v n ) relativistic 
corrections become equally important, or the ratio — proves to be of order unity. On 
the other hand, one can use the results for the short distance correction A h directly at 
the threshold region, after subtraction of the resummed Coulomb singularities. In what 
follows we discuss such a resummation. 

4. For the sake of completeness, we would like to discuss how the running of the 
QCD coupling modifies the cross section near the threshold. This issue was a subject of 
many papers. It was discussed in [^|, with the NLO accuracy and in |22]| , with the 



NNLO accuracy. Our main task here is to analyze how the perturbative expansion works 
in the scheme with the running QCD coupling and how accurate are the predictions, 
made in this scheme. We perform a numerical analysis in the position space and compare 
its results with those obtained in the fixed normalization point scheme. Our analysis 
is complementary to that made in [22 1, where a comparison of the potentials in the 
momentum and position spaces was discussed. 

Near the threshold, the ratio ^ is not small and one has to resum all 0(^(^) n ^J terms 



2jJ . For the tt production, we should also include the finite width of the t-quark, T t . 
The cross section of the process e + e~ — ► ti near the threshold normalized to the cross 
section for e + e~ — > is[] [0: 

247T 

R = N c e 2 lmG{E + T\;0,0), (14) 

s 

where G(E + iT t ; r, r') is the Green function for the Schrodinger equation: 

v 2 

(H — E — iT t )G(E + iTt] r, r') = 5(r — r ), H = — + V(r). (15) 

m 



what follows, we disregard relativistic corrections. 
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Thus, the problem of the soft corrections resummation to all orders in 24 , reduces to a 
calculation of G(E + iT t ; 0, 0) for the Schrodinger equation with a heavy quark-antiquark 
potential V(r). In the NNA approximation, the short-distance correction A h derived at 
the previous section, enters into R through the factor 1 + A h . 

It is well known, that at LO and NLO M no choice exists of the normalization 
point for the coupling constant, which would result in an appropriately similar behavior 
of the threshold cross sections calculated in the fixed point perturbative potential and 
the "running" one, respectively This means, that an account of running is actually 
important, since it allows one to resum large logarithmic corrections which arise in the 
fixed point perturbation theory. However, an account of the NNLO corrections makes two 
cross sections much closer to each other, if the high normalization point is used, fi ~ m t . 
The difference becomes an effect of C(a s (/i) 3 ) and appears to be not very large, at the 
level of 5-10%. 

Let us now discuss a scheme with the running QCD coupling. Usually a starting 
point for the analysis is a perturbative potential in the scheme with a fixed normalization 
point, in the momentum space. The "running" potential is obtained by resumming all 
logarithmic terms to the running coupling: 



2 2 



V(r)= - ^ ^ 11 + — i ^ — «i + — i ^ — J [A,y + «2 (10) 

fcyJMr'\\'\ „ . . / 5 \ 7T 2 

+ 



47T 



16/3 3 C(3) + fa (3/^1 + y + a 3 



Here r' = re 7B , is the Euler constant. The coefficients a\,a,2 are known p3|, |24| , 
while 03 is still unknown. The coupling a s (l/r') suffers from the Landau pole, appearing 
as a consequence of extrapolation into the strong coupling domain of the perturbation 
theory result . The usual way to handle this singularity is an introduction of some model 
potential, mimicking non-perturbative effects, at the large distances r > r ~ 1/Aq CD 
||, [J. However, the contribution of the "non-perturbative" region to the resulting cross 
section proves to be extremely small, less than 0.1 — 0.3%, due to the large width of the 
top quark. Therefore this cross section can be considered as model-independent. 

To find the cross section for various potentials, we solve the Schrodinger equation 
numerically. Our result for the running NNLO potential agrees with the corresponding 



result of |22[, obtained in coordinate space, when the same input parameters are chosen. 
In Fig. 2, we plot the cross section R as a function of the energy E. We have chosen 
m t = 175 GeV, a s {m z ) =0.118, Y t = 1.43 GeV. For a conservative estimate we have 
chosen a 3 = 100a 2 (recall that a 2 ~ 70ai). 

Unfortunately, the resummation of the logarithms does not place the higher order 
corrections under control. One can see from Fig. 2, that the corrections to the cross 
section are not small in the vicinity of the peak. Moreover, these corrections do not show 
a decrease in their values with increasing order. The weak convergence of the perturbative 
series is caused by the enormously large non-logarithmic coefficients, a n , entering the 
potential We see that the IS peak is shifted by about 0.5 GeV upon inclusion of 

each new order into the potential. At least partially, this effect can be due to the use of 
the pole mass. It is known that the latter suffers from the renormalon ambiguities and 
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Figure 2: R(e + e" -> tt) for the LO, NLO, NNLO, NNNLO potentials in the scheme with 
the running coupling as a function of E, GeV. We use m t =175 GeV, a s (mz) =0.118, 
T t =1.43 GeV. 



thus cannot be determined with accuracy better than Aqcd 0! P5L [T2] . The same order 
ambiguities are present in the interaction potential, V(r). As was shown in [p7|, |2^, |30|, 
the leading order O(Aqcd^) ambiguities cancel in the combination V(r) + 2m. 

It is rather natural to anticipate that a more safe way to construct an expansion is 
to use a running mass at the low normalization point fi, but the issue requires a more 
detailed study. 

Let us now comment on the accuracy of results for the cross section. Fig. 2 shows, 
that the conservative estimate of 03 suggests the deviation of the NNLO curve by about 
10% at the peak. We have also calculated the cross section with the NNLO potential in 
the scheme with the fixed normalization point and have compared it with the "running" 
NNLO result from Fig. 2. The difference is again about 10%. Our estimate of the accuracy 
confirms the conclusion of the paper . 

5. We have studied QCD radiative corrections, originating from the running of 
the coupling a s , to the total cross section of the process e + e _ — > tt near the threshold. 
We have resummed potentially large C^(/3oa s ) n ) short- distance QCD corrections to the 
production vertex ^*ti. The exact in /3oa s short-distance correction renormalizes the 
current to —9.5% with respect to the Born result and is well under control. The long 
distance corrections are accounted for by using a potential with the running QCD coupling 
and are more significant. We demonstrate that the accuracy of present predictions for the 
cross section in the scheme with the running QCD coupling is about 10%. 
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